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INFORMATION-Book Shelf 
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[1] Elementary Linear Algebra with DERIVE , R J Hill and Th A Keagy 
Chartwell - Bratt, 1995, ISBN 0 86238 403 6 

(This is a textbook with DERIVE fully integrated. DERIVE is used as a natural tool, 
370 pages) 

[2] Technology in Mathematics Teaching, L Burton and B Jaworsky, ed. 

Addison - Wesley, 1995, ISBN 3 89319 857 1 

(This is not a DERIVE- book. It addresses general issues raised by the introduction of 
technology into the teaching and learning of mathematics, 496 pages) 

Two "classics" of DERIVE literature are now available in a German translation: 

[3] Mathematik entdecken mit DERIVE - von der Algebra zur Differentialrechnung, 

Jerry Glynn, Birkhauser, 1995, ISBN 3 7643 5001 6, 1 54 Seiten 

[4] Mathematik lernen mit DERIVE, J Berry, A Watkins, E Graham, Birkhauser, 1995, 

ISBN 3 7643 5193 4, ca 360 Seiten 

[5] Introductie DERIVE 3.0, Paul Drijvers, STAM Techniek, 1995, ISBN 90 401 0640 1 
(A Dutch DERIVE book from Paul, a DERIVIAN from the first days.) 


Exchange for DERIVE Teaching 


Borse fur DERIVE Unterrichts- 

materials in the DNL 


materialien im DNL 

The wheel has not to be invented twice. 


Das Rad muss nicht zweimal erfunden werden. 


I can offer: 

Binomial Theorem, GCD & LCM, System of Coordinates (in English and in Ger- 
man as well) 


Second North American DERIVE User Group Meeting 

I have the pleasure to announce the 2 nd DERIVE User Group meeting in the 
United States of America. We will meet at the occasion of the 

ICTCM 1995, November 16-19, 1995 in Houston, Texas. 

Jerry Glynn from MathWare in Urbana, IL, has promised to help organizing 
this meeting. If you will not receive more details about the DUG meeting until 
the beginning of the conference you will find information at the booth of Math- 
Ware. I am very glad that - supported by Soft Warehouse Hawaii - 1 am able to 
be present at this event. And I hope to meet many of our American members 
and the European members, who will attend this important conference. 

I am looking forward to meeting you in Texas. 

Josef 
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Liebe DUG-Mitglieder, 

Da ich diesen Newsletter rasch fertigstellen 
mochte, wird mein Letter kilrzer als gewohnlich 
ausfallen. 

Der Sommer war ausgefullt mit DERIVE- und 
anderen Reisen. Aufierdem musste die deutsche 
Ubersetzung des Bucks "Learning Mathematics 
through DERIVE" abgeschlossen werden. 

Das herausragende Ereignis des Sommers war 
sicherlich die International DERIVE Conference 
in Honolulu. Viele der 35 Derivianer aus aller 
Welt nahmen die Konferenz zum Anlass, diesen 
wunderschonen Teil unserer Erde langer zu besu- 
chen. Und ich denke, dass keiner seinen Aufenhalt 
auf der einen oder anderen Insel bereut hat. Meine 
Frau und ich besuchten seeks Inseln und hatten 
einen unvergess lichen Urlaub. An dieser Stelle 
herzlichen Dank im Namen aller teilnehmenden 
DUG Mitglieder an Soft Warehouse Hawaii fur 
die Vorbereitung und Unterstiitzung der Konfe- 
renz, sowie fur die Gastfreundschaft und Betreu- 
ung wdhrend der Tage in Honolulu. Mahalo nui. 

In diesem Herbst wird bestimmt wieder ein Mei- 
lenstein in der Entwicklung der Computer Algeb- 
ra gesetzt. Mit dem TI-92 kommt ein " Taschen - 
rechner” auf den Markt, der ein CAS, sehr ahn- 
lich DERIVE, mit alien seinen Vorziigen imple- 
mentiert hat. Aber nicht nur das, auch CABRI- 
Geometre, ein Texteditor und eine Tabellenkalku- 
lation - die nicht nur numerisch. sondem auch 
symbolisch arbeiten kann -finden sich in diesem 
Zauberkasten, der in USA ca 200 $ kosten soil. 

Ich habe DERIVE fur WINDOWS gesehen ! 

Teresa, eine charmante und attraktive junge Da- 
me, arbeitet emsig daran, DERIVE so wie wir es 
alle schatzen, in die WINDOWS Umgebung zu 
iibertragen. Daneben ist sie eine aufierst begabte 
Hula-Tanzerin. Aloha Teresa, aus Osterreich und 
herzliche Griifie an Deine kongenialen Tanzpart- 
nerinnen. 

Barbel Barzel hat angekiindigt, dass ein ausfuhrli- 
cher Tagungsband der DERIVE Days Dilsseldorf 
herausgegeben wird. Sie werden sicher im nach- 
sten DNL daruber informiert. 

Herzliche Griifie 



Dear DUG Members, 

This letter will be shorter than usual because I 
want to finish the Newsletter in time. 

Summer was spent with DERIVE- and other 
journeys. Moreover the German translation of 
"Learning Mathematics through DERIVE" had 
to be finalized. 

Certainly the most important event of this 
summer was the International DERIVE Con- 
ference in Honolulu. On this occasion many of 
the 35 Derivians from all over the world visited 
this wonderful part of our globe. And I think 
that nobody has regreted his stay on one or 
the other island. My wife and I were on six 
islands and it was an unforgetable holiday. 

I want to express our special thanks to Soft 
Warehouse Hawaii for preparing and support- 
ing the conference and for the hospitality and 
care during our stay in Honolulu on behalf of 
all DUG members present at the conference. 
Mahalo nui. 

In this fall will be set another milestone in the 
evolution of Computer Algebra. With the TI-92 
we will have a pocket calculator with a regular 
DERIVE implemented. Additionally CABRI- 
Geometre, a text editor and a spread sheet - 
with numerical and symbolic abilities - can be 
used on this magic box. For the US a market 
price of about 200$ is announced. 

I have seen DERIVE for WINDOWS!!! 

Teresa, a really charming and attractive lady 
- is busy to transfer all that what we are ap- 
pre- dating in DERIVE into a WINDOWS en- 
viron- ment. Besides that she is a very gifted 
hula dancer. Aloha Teresa, from Austria and 
best greetings to your congenial hula part- 
ners. 


Best regards 
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The DER1VE-NE WSLE TIER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a con- 
tents of 30 pages minimum. The goals of 
the D-N-L are to enable the exchange of 
experiences made with DERIVE as well as 
to create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone: 43-(0)2275/8207 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the D-N-L. 
It must be said, though, that non-English 
articles are very welcome nonetheless. 
Your contributions will be edited but not 
assessed. By submitting articles the author 
gives his consent for reprinting it in D-N- 
L. The more contributions you will send, 
the more lively and richer in contents the 
DERIVE Newsletter will be. 


Preview: (Contributions for the next issues): 

Stability of systems of ODEs, Kozubik, SLO 

Graphic Integration, Probability Theory, Linear Programming, Bohm, A 

LOGO in DERIVE, Lechner, A 

DREIECK.MTH, Wadsack, AUS 

IMP Logo and Misguided Missiles, Sawada, HAWAII 

3D Geometry, Reichel, AUS 

Parallel- and Central Projection, Bohm, AUS 

Muller's Method to solve univariate equations a.o., Speck, NZL 

Vector and Vector Indices Sorting, Biryukov, RUS 

Algebra at A-Level, Goldstein, UK 

Tilgung fremderregter Schwingungen, Klingen, GER 

Utility for Complex Dynamic Systems, Lechner, A 

Some Improvements on the Resolution of ODEs, Fuster, E 

Notes on DERIVE 2.6 functions and limits, Speck, NZL 

and 

Setif, FRA; Vermeylen, Belgium; Lymer, FRA; Leinbach, USA, Aue, GER; 
Weth, GER, Wiesenbauer, AUT; Keunecke, GER 
and messages from the derive-news@mailbase.ac.uk 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag.Josef Bohm 
Herstellung: Selbstverlag 
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DERIVE 3.x. Bugs or not 

Preparing this issue I had three times the problem, that functions which did very well with versions 2.5x led to a 
desaster using DERIVE 3.x. (Spline-function and Bezier- functions). I found similar problems with some own 
functions, unfortunately once when I tried to demonstrate a difficult file during a presentation. During my stay 
in Honolulu I found some occasions to talk with Albert Rich about these "incompatibilities". One problem 
could be solved immediately. See the explanation in the Bezier contribution concerning the MB-function, page 
3 1 . This case will not occur very often. The other one has much more chances to happen. Back home from the 
holidays I found a letter from Albert Rich which I will cite: 

Enclosed is the DERIVE XM version 3.06 diskette I promised to send you. It fixes the bug 
caused by the user-defined function that tried to extract the 0 th element of a vector. For example, 
simplifying 0*element ( [x, y] , 0) or 0* [x] sub 0 used to hang DERIVE before a version 3.06. 

Concerning the problem with SPLINE. MTH: Often a dummy variable needs to be introduced 
in a function's definition that is not a formal argument of the function. For example, the following func- 
tion that sums the square of a vector's elements introduces the dummy variable i: 

sumsqr(v) := sum((v SUB i) A 2, i, 1, DIMENSION (v) ) 

If one function calls another and both functions introduce the same variable, the introduced 
variable's value may not be what was intended. Therefore, introduced variables in functions that call 
each other should be given unique names. 

The problem with spline. MTH was caused by using z_ as the introduced dummy variable in 
several function definitions. The revised SPLINE. MTH included on your diskette corrects the prob- 
lems by using unique variable names for introduced variables. 

Let me know as soon as you or a DUG member finds a new or old problem with DERIVE. 

Many thanks, Albert. By the way I have to mention that Albert is very interested in the DNLs as one of his most 
important sources for improving DERIVE. So go on, hunting for problems and bugs, Albert is waiting. 

Nevertheless, I cannot find it very userfriendly if functions loose their functionality by upgrading DERIVE. 
That might cause very surprised eyes during a demonstration of DERIVE's capabilities. Both changes men- 
tioned above are surely justified by mathematical exactness. But I hope for us all that those changes will be at 
least announced in a README-file. 


What you have to do now with the spline-file on page 1 8 is the following: Substitute the variable z_ in the func- 
tions SPLINE-D (e.g. by a ), SPLINE-B (by b_) and SPLINE AUX (by c_). I could have done it for you. But I 
think that you have to see this possibility of making yourself a lot of troubles and how to resolve them. (The 
diskette of the year will contain the revised version). Is obsolete with DERIVE 6. 10, the files are ok, Josef 


Wolfgang Propper. Niirnberq. Germany 

During the last weeks - testing the new TI-92 - 1 found an interesting question: 

DERIVE seems to dislike the product Va -4b .There is no trick to bring it into the form Vab. This 

would not be so bad, but if I consider the product Vl-x 2 -.E-Ly which is defined for -1 < x < 1, and 

V 1-x 

can be simplified by hand using common sense to (1+x), it becomes problematically. DERIVE refuses 
all attempts simplifying this expression. The plot is correct only for -1 < x < 1 . The plot shows -\1 + x| . 


The same will occur with expression 


U-l 


DNL: Try to substitute a = -2, b = -8, then you will find that, following the rule drawing roots before multi- 
plying expressions, there is no identity in your expressions: V- 2 • V- 8 = iyfl • z V 8 = -4 but yj (-2) • (-4) = 4. 


If you declare the variables nonnegative then you are able to proof the identity. The same idea can be applied to 
your other expression. I don’t understand, why in this case DERIVE fails to show the identity for x e [-1,1). 
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This was DERIVE 3.xx 



This is DERIVE 6.10: 

#i: -JWy = -Jx'-Jy 

#2: -Jx-Jy - 7(x-y) = -Jx-Jy - ■JC.x-y') 

#3: x :e Real (0, «) 

#4: Vx-Yy = yWy 

but: 

#5: ,/x - ,/y - 7Cx-y) = 0 



My comment of 1995 was: I believe that the plot (DERIVE 3.xx and the value tables) are correct. If you don't 
simplify then draw the roots first followed by multiplication. That is the order of performing the arithmetic op- 
erations. I hope that my explanation is right and the answer will satisfy you and the readers. Any other com- 
ments are highly appreciated. My question about the missing result 0 in expressions #9 (DERIVE 3) and #10 
(DERIVE 6) might be answered in this way, too. 

How are the answers of the Tl-devices? Voyage 200 and TI-NspireCAS behave pretty the same: 




|.x>0 and y> 0 



It performs the simplification of the product of the roots. 
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The next screenshots show how the other expressions are treated by the TIs: 



Darlene Boharsik, Marietta, Georgia. USA 

Thank you for the timely manner in which you forwarded to me all 1 7 issues of the published DNL. 

I am enjoying every copy. I look forward to all future issues. It is quite interesting I find, reading the 
forum from users around the world. 

I am working towards my second career as a math teacher and discovered DERIVE 2 years ago. 

I feel that there is a need for the US schools to implement a program such as DERIVE as an addi- 
tional tool in the math classroom. I hope to help bring this knowledge to our math departments. Our 
math department at the Southern College of Technology at this time is a bit slow on change. 

Please keep up the good work as I will be spreading the news about DERIVE's use in the class- 
room. 


Andreas Zwicknagel, Regensburg, Germany 

Lieber Josef, dear members of the DUG! 

Schon vom ersten Durchblattem der DNLs bin ich tief beeindruckt. Ich habe gar nicht gewusst Oder 
geahnt, zu was DERIVE fahig ist. Zudem bewundere ich das Engagement des Editors. Wir beide ha- 
ben zwar denselben "Nebenberuf", aber die Zeit zu einem tiefer schurfenden Einstieg in die Pro- 
grammierfahigkeiten anhand der DNLs wird mir wohl erst in den nachsten Sommerferien bleiben ! 

Mit freundlichen Gruften, A.Z. 


DNL: I hope you had some time to work with DERIVE during the holidays. For you and all our members I add 
one plot which I have received from a young and very enthusiastic Belgian DERIVE - math teacher. I met him at 
Turnhout. where I had a workshop about graphics with DERIVE. He told that his students -s econdary school - 
use DERIVE to design phantasy buildings. He promised to send some prints. Dear Luk Van den Broeck, thank 
you very much for not only promising, but promptly sending. Congratulations to your students. You can find 
other spread in this issue. 



More contributions to the User Forum on page 45 
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Didactics of CAS 
Conic Sections Escape R 2 

Dr. Karl Fuchs, ACDCA, Salzburg, Vienna 


0 Didactical Introduction 

In the 7 th class of the Austrian Grammar school we have to teach the conic sections: ellipse, hyperbola 
and parabola. Although we have the methods of analytical geometry to create the equations of conic 
sections in R 3 they have no chance to escape the R 2 in the moment. We are discussing focus-characte- 
ristics of conic sections: 

XF x + XF 2 = 2a; 2 a > F X F 2 and e 2 = a 1 - b 2 for the ellipse, 

| XF x + XF 2 1 = 2a; 2 a < F X F 2 and e 1 = a 2 + b 2 for the hyperbola, 

XI = XF for the parabola, 


derive the equations: 

b 2 x 2 +a 2 x 2 =a 2 b\ b 2 x 2 -a 2 x 2 = a 2 b 2 and y 2 =2px, 

intersect the curves and even intersect them with straight lines. But although we are doing such a lot 
of things with the curves the students ask themselves: ’’Why are these curves called conic sections?” 
Now the teacher has the possibility to refer the student to a picture in the textbook that shows a cone 
intersected by different planes generating the discussed curves and the teacher can continue telling 
them that this is the reason why these curves are called conic curves. But most of our students - 
especially nowadays - aren't interested in telling, they are interested in doing things - even mathe- 
matical things. 

In 1983 Univ.Prof. Dr. Fritz Schweiger, Professor for the Didactics of Mathematics at the University 
of Salzburg wrote a didactical article dealing with the problem of teaching conic sections in Grammar 
school. In an appendix he suggested to escape to R 3 to make these curves interesting for students with 
high ability in mathematics. At that time we were able to do the analytical geometry in mathematics, 
the illustration of the processes and results we had to pass on to the teacher in geometry with regard to 
construction which is an independent subject called ’Darstellende Geometrie' in Austrian Grammar 
schools. This was really a problem. 

I think nowadays with DERIVE we have an instrument to do these illustrations, that we better call 

visualizations now because of the dynamic aspect of drawing on the computer screen. ” DERIVE isn’t 

good at plotting 3D-objects” is an argument I often hear from several different users of DERIVE - 
especially from teachers. Certainly on one hand we don’t want DERIVE to go the way of the huge 
CAS-systems like MAPLE or MATHEMATICA in graphics, which are too big for the school but there 
is a lot we can do with DERIVE - even with 3D-objects. 
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1 The Problem 

The two objects of our investigations are 
cone k with a right angle at the top 


( 1 ) 

( 2 ) 


: z 2 and 


k\x 2 + y 2 

an intersecting plane e containing point 
i\0,0,zi), parallel to the x-axis with an 

angle 9 of inclination 


e:< 


x-u 

y = v-cosO 
z = Zj + v • sin 0 


a = (1, 0, 0) _L b = (0, cos 9, sin 9) 



The figure shows the two objects in a front view and a top view (horizontal and frontal normal projec- 
tion). 

Let's find the general form of the curve in plane e intersecting the cone k and the plane e. 
kne : 

u 2 + v 2 cos 2 0 = z 2 + 2 Zj v sin 0 + v 2 sin 2 0 
u 2 + v 2 (cos 2 0 - sin 2 0) - 2 z l v sin 0 = z 2 
u 2 +v 2 cos26^-2z 1 vsin6^ = z 1 2 

The result is the equation of a curve in plane e. I decided to perform this sequence by hand to show the 
students that there is no miracle in manipulating the equation (elementary algebraic transformation and 
the trig identity cos 2 0 - sin 2 0 = cos 29) in order to obtain the general form of the curve in e. 


1 The DERIVE Session 

Now we are going to visualize the curves for different values of 0 (using functions of utility file 
GraphicsFunctions . mth). We choose z\ = 0.5 for the following examples. 

1 st step: Defining the curve and deriving the functional form u =fly , 9) for plotting. 

#1: L0AD(C: \Programme\TI Education\Deri ve 6\Math\GraphicsFunctions.mth) 

2 2 2 
#2: u + v • C0S(2 • 0) - v-SIN(0) = 0.5 

2 2 2 
#3: S0LVE(u + v .COS(2-0) - v-SIN(0) = 0.5 , u) 

2 

V(- 4* v • C0S(2 • 0) + 4 • v* SIN(0) + 1) 

#4: u = - v u = 

2 

2 

VC- 4* v • C0S(2 • 0) + 4-V-SINC0) + 1) 


2 





p 
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Our solution vector will have the form \u, v cos (0), Zi + v sin(#)]. 


2 nd step: Finding the curve for 6 = 0. 

2.1 Substitute 0 for 6 in #4 and simplify: 


#5: 


# 6 : 

2.2 

#7: 

#8: 


2 

VC- 4-v -COSC2-0) + 4-v-SINCO) + 1) 

u = - v u = 

2 

2 

V(_ 4 . v .COSC2.0) + 4.V.SINC0) + 1) 


2 

2 2 
V(1 - 4-v ) V(1 - 4-V ) 

u = - V u = 


2 


2 


Create vectors #1 and # 2 as solution vectors: 


kl := 


k2 := 


2 

7(1 - 4-v ) 

, v-COS(O) , 0.5 + 

2 


2 

7(1 - 4* v ) 

, v-COS(O) , 0.5 

2 


v-SIN(O) 


+ v-SIN(O) 


3 rd step Call the axes (contained in the preloaded utility file) and apply the isometric projection 
on k\ and k 2 \ 

#9: axes 

#10: [ISOMETRIC(kl) , IS0METRIC(k2) ] 


4 th step: Plot #9 with the parameter between -1 and +1. 

Plot #10 with the parameter -0.5 < v < +0.5 (Take care for the domain - expr #7 & #8! !) 

It might be a good question for students, why -n < v < +7i does not result in a complete plot! 



We have the isometric picture of a circle together with the axes. 
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We can generate the other curves in a very similar way: 
(I redefine the general form of the intersection curves.) 


#12: k_l 


VC- 4 ■ v ■ COS( 2 ■ 0) + 4 ■ v ■ 5IN( 0) +1) 1 

, V-COSC0), V-SINC0) + 

2 2 


#13: k_2 


VC- 4 ■ v -COSC2-0) + 4-V-SINC0) +1) 1 

, V.C0SC0), V-SINC0) + 

2 2 


71 


We obtain the intersection as an ellipse for # = — : 


#14: SUBST 


[ISQMETRICCk_l) , ISOMETRICC k_2) ] , 0, 

6 


#15 : 


2 2 
VC- 2-v + 2 ■ v + 1) VS-y VC- 2-v + 2 ■ v + 1) 

+ + v 


2 2 

VS-v VC- 2-v + 2 ■ v + 1) VC- 2-v + 2 ■ v + 1) 


1 V^ 

L 2 4 


i. - ill . -L 

L 2 4 J 2 


2 1 Vs Vs 1 

#16: SQLVE(- 2-v + 2-v + 1 > 0, v) = < v < + 

2 2 2 2 


#17: axes 



The solutions of inequality #16 are needed for entering the boundaries of the parameter v when plot- 
ting #15. It might be a good question for students why plots generated with larger ranges for v, 
eg -5 < v < +5 show gaps in the plots! ! (Try it! The square root must be entered as sqrt(..).) 


71 


We obtain the intersection as a parabola for 6 = — . Why? 


#18: SUBST 


[ ISQMETRICC k_l) , ISQMETRICC k_2) ] , 0, 

4 


#19: 


VC2-V2-Y + 1) V2-v VC2-V2-V + V2-v 1 

+ + + 

2 2 4 4 2 

V2-v VC^'V^-v + 1) VC^-V^'Y + 1) V^-y 1 

- - + + 

2 2 4 4 2 


#20: SQLVEC2-V2-Y + 1 > 0, y) = v > - 


U 


4 


plO 
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The plot isn't very attracting und you cannot really see the parabola. To obtain a better visualization 
the object was rotated applying ROTATE_Z ( 71 / 4 ) - which results in a 90°-tum. (Take care for trans- 
posing the k_l- and k_2 -vectors before multiplying with the rotation matrix!) 


#21: SUBST 


ISOMETRIC 


ROTATE_Z 


n 

k 4 


[ k_l ] 1 , ISOMETRIC 


ROTATE_Z 


n 

k 4 


[ k_Z ] 


, 0. 


# 22 : 


+ 1) ^2-y 1 ' 

v , + + 

4 2 2 . 

<]2'4(_2'42-v + 1 ) *J2 ■ v 1 

v , - + + 

4 2 2 



#23 : 


SUBST 


[ ISOMETRIC( k_l) , ISOMETRIC( k_2) ] , 0, 



#24: 


2 

-/( 2 ■ y + 2-Jl‘V + 1 ) 


2 



2 

V(2-v + 2-V^-v + 1) 


4 




2 

v -/(Z.v + 2 -^B-v + 1) 

2 2 


2 

V(2-v + 2 -^ 3 -v + 1) 


4 


+ v 




2 

#25: S0LVE(2-v + 2-^3-v + 1 > 0, v) 


V3 

v < - 


1 1 v' 3 ' 

V Y > - 


2 


2 


2 


2 



n 


One case is left. Intersecting the cone with a plane through its vertex (with 6 = — ). The plane is verti- 
cal now. The intersection curve is a pair of straight lines (real or complex). 


#26: SUBST 


[ ISOMETRIC( k_l) , ISGMETRIC( k_2) ] , 0, 


#27: 


2 ■ v + 1 


|2-v + 1| 


| 2 ■ v + 1 | 1 

+ v + 

4 2 

| 2 ■ v + 1 | 1 

+ v + 

4 2 



Here again the plot is not very informative. 

I took another plane through the vertex (z\ = 0 and 0 = 3 tt/8): 

2 2 

#28: u + v -CO5(2-0) = 0 

2 2 

#29: SOLVE(u + v -005(2-0) = 0 , u) 

#50: u = - | v | -VC- CO5(2-0)) v u = | v | ■ V(- CO5(2-0)) 

#51: col := [- |v|-V(- COS(2-0)), v-COS(0), v-SIN(0)] 

#52: co 2 := [|v|-^(- CO5(2-0)), v-CO5(0), v-5IN(0)] 

f 5-tt ^ 


#55: SUBST 


[ISOMETRIC(col) , ISOMETRIC( co2) ] , 0 


8 
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# 34 : 


3/4 3/4 

2 • | v | v.J(2 - J2) 2 ■ | v | v.J(10 - J2) 

+ + 

2 2 4 4 


v^2 - </2) 


3/4 3/4 

2 ■ | v | vVClO - ^2) 2 ■ |v | 


2 


4 


4 



References 

Schweiger Fritz , Analytische Geometrie oder sind die Kegelschnitte noch zu retten? 
Mathematik im Unterricht, Nr. 7, 27 - 29, 1983 


3D plots of Karl's Conic Sections 

Josef Bohm, Wurmla, Austria 

In July I was invited by Jan Vermeylen to give a presentation and to hold a workshop at the 
Conference of the VVWL (Vlaamse Vereniging Wiskunde Leraaren = Association of Flemish 
Mathematics Teachers) at Turnhout, Belgium. Jan gave a workshop “Conic Sections” and 
asked to produce 3D-pictures of the curves. I tried and I add the plots to Karl's examples 
from the contribution above using my results. 

These are the results with Karl's special assumptions. In the following you can find a more 
generalized solution of the problem. I must admit that it took some time to find appropriate 
values for the eye-coordinates in order to obtain nice looking pictures. 

The Circle 



y 
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The Ellipse 



The Parabola 
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I define the plane E by its intercepts on the axes and the cone: 

c_-(b_-x + a_-(y - b_)) 

#1: E(a, b, c) := lim lim lim - 

c_->c b_->b a_->a 


a_*b_ 


2 2 

#2: con := 7(x + y ) 


#3: C0SECC6, zl) := MIN 

#4: COSECCO, 0.5) 


#5: 


2 2 

2-7(x + y ) + 1 




r 

zl 


con , 

, E 

00 , - 

, zl 

V 



TAN(0) 

JJ 



2 

2 



2-7(x 

+ y ) - 1 



4 4 

The settings for the plot box (circle) are: 


#6: COSEC 


#7 : COSEC 


TT 

, 0.5 

6 


7T 


, 0.5 


*1 

Minimum 

Maximum 

Scale 


|075 

|0.75 

y:l 




f075 

|075 

|0.75 

=:l 




|3i 

h 

0.55 






#8: COSEC 


#9: COSEC 


TT 

, 0.5 

3 


TT 


, 0.5 


#10: COSECC1.57, 0.5) 


#11: MIN 


con, 0*x + TAN 


3 • 7T 


8 


n n 

In times of Derive 6 we introduce a slider bar for parameter t with — <t< — . 

2 2 


#12: C0SEC(t, 0.5) 
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The more generalized version: Find the intersection curve between the coneO(h,r) and the 
plane E(a,b,c) with its intercepts on the axes a, b and c. 

2 2 

r«h - 7( x + y ) -h 

#1: CONEO(h, r) := 

r 


c_*(b_*x + a_-(y - b_)) 

#2: E(a, b, c) := lim lim lim - 

c_->c b_->b a_^a a_*b_ 

#3: C0NESEC(h, r, a, b, c) := MAX(MIN(CONEO(h , r) , E(a, b, c)), 0) 
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#8: C0NESEC(6, 3, 2.5, 6) 


#9: C0NESEC(6, 3, a, <», c) 


Cubic Splines Part 2 & 3 

Otto Reichel, St. Polten & Leo Klingen, Bonn 

In the second part of the contribution dealing with Cubic Splines you can see how Otto Reichel solved 
the task. There are some similarities with my approach. We talked about and very soon we found the 
reason for this fact: we used the same books about computer graphics. Nevertheless I think that it 
could be interesting to compare and to find the differences. I am sure that you can follow the ‘pro- 
gram’ without further explanations if you use the theoretic basics about Cubic Splines (1). Otto used in 
his files the same nodes as I did. Then he added a more extended example. I add the plot. 


# 1 : 

#2: 

#3: 

#4: 

#5: 

#6: 


Cubic Splines by 0. Reichel; RSPLI.MTH 

H(a) := VECTORCa - a , k, DIM(a) - 1) 

k + 1,1 k,l 

Hl(a) APPEND([0] , H(a)) 

F(a) :z VECT0R(2.((H(a)) + (H(a)) ), k, 2, DIM(H(a))) 

k k - 1 

Fl(a) APPEND([0] , F(a)) 


G(a) := APPEND 


[0], VECTOR 

r 

3* 

r a - a 

k + 1,2 k, 2 

a - a > 

k, 2 k - 1,2 

(H(a)) 

(H(a)) 

V 

V 

k 

k - 1 J 


k, 


2 , 


#7: 


DIM(a) - 1 

J j 

MATl(a) := VECT0R(VECT0R(IF(i • j = 1, 
0, IF( | i - j| =1, (Hl(a)) 


MAX(i , 


1 , 

J 

j) 


IF(i = j, CFl(a)) , IF(i = 1 v j = 1, 
i 

0)))), i, DIM(a) - 1), j, DIM(a) - 1) 


-1 

#8: CK(a) MATl(a) -G(a) 
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You can immediately plot expression #20; don't forget to activate Options > Approximate before plot- 
ting in the 2D-Plot Window. 



#22: example2 := 


0 0 ' 
1 1 
2 3 

4 1 

5 0 

6 2 

7 5 

8 1 

9 3 


10 0 


Perform the same procedure as above, leading to KUBSPLI(example2): 



You have to define the result of CK(example2) as c:= and then plot 
KUB SPLI(example2) . 


Leo Klingen's Splines 

Leo Klingen wrote the most compact DERIVE-program fort he Cubic Splines. As a special 
you will find a function for parametric splines, too. Additionally there is another function to 
approximate any given function by cubic splines. 

(Even in 1995 some problems occurred with this file changing from one DERIVE-version to 
the other. I adapted the file for DERIVE 6.10. Josef) 


#1: Notation := Decimal 


rSPLINE_AL(h, z) := h , SPLINE_Bl(h, z) := 2- (h + h )1 
#2:L z - 1 z-lzj 

5PLINE_Cl(h , z) := h 
#3: z 


SPLINE_KL(a, h, z) := 3- 


a - a 
z z - 1 


h 

z - 1 j 


a_ + 1 


#5: SPLINE_D(c , h, n) := VECTOR 


3-h 


i a_, n - 1 


D-N-L#19 


Leo Klingen: Cubic Splines 3 


pl9 


#6: SPLINE_B(a, c, h, n) := VECTOR 


a -a h ■ (2- c + c ) 

b_ + 1 b_ b_ b_ b_ + 1 


R (( a , 


#7: SPLINE_AUX(m , a, c, n, h, x) := VECTOR] | a + (SPLINE_B(a, c, h, n)) .(x - m ) + c - (x - m ) + (SPLINE_D(c , h, n)) - (x - 

c_ c_,l c_ c_,l c_ 


nn ) - xCm i x , m ), c_, n - 1 

c_,l J c_,l c_ + 1,1 

SPLINE_AUX2(m , n, x) :z SPLINE_AUX(rn , m 1 , ELEMENT(ROW_REDUCE(APPEND(APPEND( [VECTOR (IF (s_ =s 1, 1, 0), s_, n) ] , VECTOR (VECTOR (I F C 
#8: 2 

-1, SPLINE_Al(APPEND(VECTOR(r[i’ - m 1 , i_, n - 1) , [0]), z_) , IF(s_ = z_, SPLINE_Bl(APPEND(VECTOR(m’ - rn ' 

l,i_ + l 1, T_ 1, T_ + 1 1, T_ 

i'_, n - 1), [0]), z_), IP<SL = z_ + 1, SPUN E_C1 ( AP PEND ( VECTO R (m 1 - rn’ , i_, n - 1), [0]) , z_) , 0))), s_, n) , z_, 2, n - 

1, f_ + 1 1, T_ 

1)), [VECTOR (IF (s_ = n , 1, 0), s_, n)]), APPEND([[0]] , APPEND^VECTORtSPLINE.DUm 1 , APP END (VECTOR (rn ’ - m 1 , i_, n - 1), 

L 2 l,i_+l 1 , T_ 

[0]), z_) , z_, 2, n -1)1’, [[0]])))’, n +1), n, APPEND(VECTOR(m’ - rn’ , i_, n - 1) , [0]), x) 

J l,i_ + l 1, T_ 

#9: SPLINE(m, x) := SPLINE_AUX2 (m , DIMENSION (m) , x) 

SPLINE_INTERPOLATION(f , x, a) := SPLINE (VECTOR ([a , lira f' 1 , k_, 1, DIMENSION (a) ) , x) 

#10: L k_ x -*aj, k_ J 

-4 -3 
-2 2 
0 1 

2 5 

3 4 

5 7 

[-432 -4 3 2 

#12: Ll. 31027-10 - (1025-x + 12300- x + 35560- x + 22488) - SIGN (x + 2) - 1.31027-10 -(1025-x + 12300- x + 35560- x + 22488) - SIGN (x + 4), 1.31 

#12 expanded: 

[32 3 

#13: L - 0. 134302675 -x -SIGN(x + 4) - 1. 6116321- x -SIGN(x + 4) - 4 .65932012-x- SIGN(x + 4) - 2 .946535176- SIGN (x + 4) + 0. 134302675- x -SIGN(x + 


#11: SPLINE 



#14: 


#15: 


|^SIN(x), x, | — jt , - , 0, , rrjj 


SP L I N E_I NTE R PO L ATIO N 


[ -10 6 2 7 7 -10 6 2 7 

Ll. 25926-10 -(4.53299-10 -x + 2.84816-10 -x + 1.11847-10 )-|ll3-x + 355 | - 1.25926-10 ■ (113 -x + 355) ■ (4 . 53299 ■ 10 -x + 2.84816-10 ■: 

r ^ 2 

L - 0. 06450280675 -x -SIGN(226-x + 355) - 0. 6079244187 -x -SIGN(226-x + 355) - 1.432387988-x- SIGN(226-x + 355) - 0.4999978089- SIGN (226- x + 



# 16 : 
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#17: 5 ex am pie- points for parametric plotting 


#18: 


#19 : 


# 20 : 


'10' 

1 2 
P := 3 3 

4 1 
.32. 
vCP-, O : = 

Prog 

v_ := [0] 

i := 2 

Loop 

If i > DIMCpJ 
RETURN v 

v := APPEND(v_, [FIRST(REVERSE(v_)) + ^/CCP-tul - P-tO' - 1)U) A 2 + (p_4,u2 - p_^(i' - l)j,2) A 2)]) 

i :+ 1 

SPLINE_PARAM(p_, t_) := VECTQRjT (SPLINE (|V(p_) , VECT0R(p , i, DIMENSION (p_) ) j’ , t_)) , (SPLINE (|U(p_) , VECT0R(p , 


DIMENSION (pj) 


#21: SPLINE_PARAM(p) 


# 22 : 


3 

- 0 . 0292300B489 ■ t_ -SIGN(t_ - 2) + 0 

3 2 

0. 03575783648- t_ ■ SIGN (416020- t_ - 1762289) - 0. 3899273482- t_ ■ SIGN (416 020- t_ - 1762289) + 0.8967749163- t 
3 2 

- 0. 02279711254- 1_ ■ SIGN (416020- t_ - 1762289) + 0. 22B21982B4- t_ ■ SIGN (416020- t_ - 1762289) - 0.1S90626526- t_- SIGN (416020- t_ - 1762289) 
3 2 

- 0 . 0B124B963B2 ■ t_ ■ SIGN (363401- 1_ - 2881680) + 1. 212430738- 1_ - SIGN (363401- 1_ - 2881680) - 9 .103607236- 1_- SIGN (365401- 1_ - 288168 


#23: Before parametric plotting display v(p) ; then take for t_ the interval 1 0 < t_ < last element of v(p) (here: 7.88) 

#24: v (p) 

#25: [0, 2, 4.236067977, 6.472135955, 7.886349517] 



SPLINE (matrix, variable) returns the cubic splines for the points given in matrix with 
variable variable. 

SPLINE_INTERPOLATION (function, variable, list of x-values) produces the 
cubic spline approximation of function with variable. 

SPLINE_PARAM (matrix, variable ) returns the parametric splines. The domain of the variable 
for plotting is from 0 to the last number which appears in v (matrix) . 


Mr Klingen won the competition. I don't hesitate to award him the virtual 

Golden Spline. 
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Jan Vermeylen from Belgium created this nice cat using a polygon consisting of 71 points. 
First I plot the points in discrete and then in connected mode. 



Then I used Leo Klingens PARAMSPLI to give the cat a smoother contour. I couldn't do it in 
one step - this is too much for DERIVE 's memory. I had to take the cat to pieces: 


catl cat 


#16: 

[i 

20] 

#17: 

FIRST(REVERSE (v (catl) ) ) 

#1S : 

19 . 217 


#19 : 

SPLINE_PARAM(catl l 

O 

#20: 

cat2 := cat 

[20, ... 

, 40] 

#21: 

FIRST(REVERSE (v (cat2) ) ) 

#22: 

20.24S 



You can distinguish see the four pieces by the col- 
ours. 

The right cat was filled with colour applying the 
conversion to a picture object (Edit Menu) followed 
by little work with MS Paintbrush. 


and so on. cat is the matrix of the points. 
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Ein Zugang zu den Bezierkurven mit DERIVE 
An Approach to the Bezier Curves with DERIVE 

Gunter Scheu, Pfinztal, Germany 

1. Einleitung 

Mit der Verbreitung der frei skalierbaren True- Type-Schriften wurden auch die Bezierkurven allge- 
mein bekannt, da sie auch zu deren Darstellung benutzt werden. In den sechziger Jahren haben Paul 
de Casteljau bei Citroen und Pierre Etienne Bezier bei Renault unabhangig voneinander dieses Kon- 
zept der Dar- stellung von Kurven und Flachen entwickelt. Sie haben sich als Konstrukteure in der 
Autoindustrie mit dem Problem beschaftigt, Autos am Zeichenbrett zu konstruieren, Automodelle zu 
entwerfen und Autoteile mit NC-Werkzeugmaschinen zu produzieren (ahnliche Probleme gab es 
wohl auch in der Flugzeugindustrie). Diese Produktionsschritte liebcn sich dadurch enorm vereinfa- 
chen, dass es gelang, die Kurven und Flachen einfach zu beschreiben, zu konstruieren und mit weni- 
gen Punkten abzuspeichern. 

Mit der Methode von de Casteljau und Bezier konnen beliebige Kurven und Flachen mit Hilfe der 
Bezierpolynome dargestellt werden. Diese Darstellung ist eine Approximation der gegebenen Kurve 
bzw . Flache, da diese i.a. nicht analytisch darstellbar sein miissen. Die Polynome tragen den Namen 
von Bezier, der im Gegensatz zu de Casteljau seine Ergebnisse veroffentlichen konnte (Citroen be- 
hielt die Ergebnisse von de Casteljau geheim, sie wurden erst ab 1975 bekannt). 

Wir wollen die wichtigsten Eigenschaften der Bezierpolynome aufzeigen und damit Kurven als Funk- 
tion eines Parameters entwerfen und zeichnen. 


Together with the spreading of the free scaleable True Type fonts the Bezier curves became 
familiar because they are used for presentation of these fonts. In the sixties Paul de Caste- 
jau at Citroen and Pierre Etienne Bezier at Renault developed independently this concept of 
representing curves and surfaces. As designers and constructors working in car-industry 
they dealt with the problem to design cars on the drawing board, to construct their models 
and to produce the parts of the cars using NC-machines. These production steps could be 
simplified in a high degree because it was now possible to describe curves and surfaces in 
an easy way and to save them using only a few points. 

Applying the method of de Casteljau and Bezier any arbitrary function can be described by 
Bezier polynomials. This is an approximation of the curves or surfaces wanted which have 
not to be analytically describable. The polynomials are called Bezier polynomials because 
Bezier published his results while de Casteljau was forced by Citroen to keep his work se- 
cret until 1975. 

We want to show the most important properties of the Bezier polynomials and then design 
and plot curves as functions of one parameter. We use Bernstein polynomials which have 
only one maximum value in the given interval and are passing the endpoints of the interval. 
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2. Bezierpolynome 

Bezierpolynome werden mit Hilfe der Bernsteinpolynome dargestellt, die im gegebenen Intervall nur 
ein Maximum haben und durch die Endpunkte des Intervalls gehen. 


Die Bernsteinpolynome lauten mit den Binomialkoeffizienten 




nl 








i\(n- i ) ! 

(1 t g [0,1], i-0, (the Bernstein polynomials) 


Es sind Polynome vom Grad n mit den Eigenschaften ( Polynomials of degree n with properties) 


a) = Z 




(l - t) n ~' f - ((l - 0 + 0" - 1 


b) B n i ( t ) hat eine /-fache Nullstelle fur t = 0, (a zero of order i for t = 0) 

c) Z?" ; (? ) hat eine «-/- fachc Nullstelle fur t = 1, (a zero of order n-1 for t = 1) 

d) B". ( t ) hat an der Stelle t = — sein einziges Maximum, ( sole maximum at f = — ) 

« n 


e) B n j (t) = B" n i (1 - 1), d.h. sie sind symmetrisch bez. t und t - 1. ( symmetric wrt t and t- 1) 


Wir veranschaulichen uns grafisch einige Bernsteinpolynome. 

Die Beispiele sind in DERIVE-Notation 

Visualization of some Bernstein polynomials (in DERIVE notation) 


n-i 


#1: 

B(i , 

n, t) 

:= COMBCn, 

O- 

Cl - 

t) 

#2: 

VECTOR (B (i 

, 4, t), i 

, o, 

4) 


#3: 

VECTOR (B (i 

, 5 , t) , i 

, o, 

5) 


#4: 

FCt, 

n) := 

VECTOR C[t, 

BCi 

i n i 

t)], i 

#B : 

F C t , 

4) 





#6: 

FCt, 

5) 







i 

t 


0, n) 



Die 5 Bernsteinpolynome vom Grad 4 und die 6 Bernsteinpolynome vom Grad 5 
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Die Bernsteinpolynome geniigen der Rekursion: 

A recursion formula for creating the Bernstein polynomials: 

B" i (t) = (\-t)B"-\(t) + t-B n -\_,(t) mit S° 0 (0 = 1, B n j(t) = 0 fur j > n, *e[0,l]. 

Die Bernsteinpolynome sind wegen der Eigenschaften b) und c) linear unabhangig. 

The Bernstein polynomials are linear independent because of their properties b) and c). 

Das Bezierpolynom vom Grad n ist eine Linearkombination der Bernsteinpolynome vom 
Grad n. 

The Bezier-polynomial of degree n is a linear combination of Bernstein polynomi- 
als of degree n. 

p i ,(t) = f j h r B;‘(t), ne N, fe[0,l] 

z=0 

Durch die Bezierkoeffizienten b if i = 0, . . n ist das Bezierpolynom und damit die Bezierkurve ein- 
deutig festgelegt. 

Die Bezierpolynome haben die Eigenschaften ( Characteristics of the Bezier polynomials ): 

A) In den Endpunkten des Intervalls gelten p n ( 0) = b 0 und p n ( 1) = b n . 


B) 


Da das Bernsteinpolynom 5"- bei t 


— sein Maximum hat, sind die Punkte 
n 



i = 0, . . n von besonderer Bedeutung. Das Bezierpolynom geht zwar nur durch die Endpunkte 
(0, bo) und (1, by) aber es approximiert das Bezierpolygon durch die anderen Punkte. Die 
Punkte des Bezierpolygons liegen aquidistant im Intervall [0,1]. 


The points (i/n, bj) are of special importance because they are the maximum points of 
the Bernstein polynomials. The Bezier polynomial passes only the endpoints (0, b 0 ) 
and ( 1, b n ) but it needs the other points for approximating the Bezier polygon. The 
points of this Bezier polygon are lying equidistant in [0,1]. 


C) Fur die Ableitungen in den Endpunkten des Intervalls gelten: 

p\ 0) = n(b x -b 0 ) und p\ 1) = n(b n - b n _ x ). Die beiden ersten und die beiden letzten Koeffizien- 
ten bestimmen die Steigung der Bezierkurve in den Endpunkten des Intervalls. Analoge Glei- 
chungen gelten auch fur die hoheren Ableitungen. 


Both first and last coefficients define the slope of the Bezier curve in the endpoints of 
the interval. Analogous equations are valid for higher derivatives. 


Bezier hat nicht die Bernsteinpolynome as Basis benutzt. Er stellte seine Kurven als Linearkombina- 
tionen von Funktionen dar, die die gleiche Rekursion erfiillen wie die Bernsteinpolynome, wobei sich 
nur die Anfangsbedingungen unterscheiden. R. Forrest entdeckte, dass die Bezier-Form einer Bezier- 
kurve mit Hilfe der Bernsteinpolynome dargestellt werden kann. 


Bezier didn 't use Bernstein polynomials as a base. He presented his curves as linear com- 
binations of functions, which meet the same recursion as the Bernstein polynomials but dif- 
fer in the initial conditions. R. Forrest found out that the Bezier-form of a Bezier curve can be 
represented by means of Bernstein polynomials. 
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3. Zeichnen von Kurven mit Bezierpolynomen 

Wir betrachten nur Kurven in der Ebene, d.h. Kurven, die durch zwei Komponenten als Funktion 
eines Parameters dargestellt werden. Raumkurven werden mit einer Komponente mehr dargestellt. 

We will consider curves in the plane only, i.e. curves which can be represented by two com- 
ponents as a function of one parameter. Space curves need one more component. 


Pn{t) = Y J b i B n iif\ neN,te[ 0,1]; p n :[0,1] -> R 2 , b, e R 2 . 

i=0 

Wir versuchen die Sinusfunktion im Intervall [0, 2n] durch Bezierkurven zu approximieren. 
We try to approximate the sine function in [0, 2n ] by Bezier curves, 
p ist die verwendete Punktematrix. (p ist he matrix of the given points ). 


FB(p, t) := VECT0R(p -B(W - 1, DIMENSION(p) - 1, t) , i, DIMENSION(p)) 
#9: i 

DIMENSION (p) 

#10: GB( P| t) := I (FB( P| t)) 

i'=l i 


Die Funktion GB erzeugt die Bezierkurve fiir t e [0,1]. Wir wahlen zuerst die Punktmatrix si - fur 
sin(x) - und dann s2 - fiir 2 sin(x). Mit GB ( s 1 ) bzw GB ( s2 ) erhalten wir die entsprechenden Be- 
zierkurven. 


#9: SIN(x) 


#10: si := TABLE 


SINCx), x, 0, 2 -7T, 


#11: FB(sl) 


3 3 

2-n , -t-(l - t) 4-t-(l - t) 

2 2 

6-rr- 1 ■ (t - 1) 0 

3 3 

6-rr-t ■ (1 - t) 4-t ■ (t - 1) 


2 -^^ 0 
#12: GB(sl) = [2 -tt- t , 4-t-(t - l)-(2-t - 1)] 


#13: 


s2 := TABLE 


2- SINCx), x, 0, 2 -tt, 





#14: GB(s2) = [2 - jt- t, - l)-(2-t - 1) ] 
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GB gives the Bezier curve for 0 <t < It. We choose point matrices si and s2 and obtain the 
two curves calling functions GB(sl) and GB(s2). 

Nun wollen wir versuchen, die Konturen von zB Vasen Buchstaben oder Ziffem mit Hilfe zweier 
Bezierkurven zu erreichen. Dazu wahlen wir die beiden Punktmatrizen q\ und q2 und zeichnen die 
Kurven wieder mit 0 < t < 1 . 

Now we will design the contours of eg vases, letters or numbers by two Bezier curves, 

(ql and q3 make a "3") and plot with 0 <t <1. 


#15: 


" 

' 1 

0 ' 



" 





' 1 

8 " 



4 

2 





ql := 



1 q2 := 

4 

10 



4 

6 









. 1 

10 . 



. 1 

8 . 





#16: [GB(ql, x) , GB(q2 , x)] 


2 2 

- 9-x + 9-x + 1 - 2-x- (2-x - 3-x - 3) 

2 2 

- 6-x + 6-x + 1 - 2-(x -2-x- 4) 



(Wenn man die Punkte unterschiedlich "gewichtet" (= mehrfach in die Matrix aufnimmt) andert sich 
der Kurvenverlauf deutlich. Josef) 

(If one gives the points different "weights" (= by taking them several times into the matrix) 
the run of the curve changes significantly. Josef) 

#17: GB(APPEND(ql, q2) , x) 


#18 : GB 


r- 

" 1 

0 " 



4 

2 



4 

6 

. x 


1 

8 



4 

10 


X 

. 1 

10 . 

J 


#21: 


[210.) 


#15: GB 



" 1 

0 " 



4 

2 



4 

6 



1 

8 

■ x 


1 

8 



1 

S 



4 

10 



. 1 

10 . 

J 


#20: GB 



" 1 

0 " 

N 


4 

2 



4 

6 



1 

S 



1 

S 



1 

S 

1 X 


1 

S 



1 

S 
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10 9876542 

- 1050 • x + 2835 • x - 5040-x + 5670-x - 3780-x + 1260-x - 135-x + 


9 8 7 6 5 4 3 

30 • x + 1, 2 • x • ( 17 • x + 50 • x - 315-x + 480-x - 210-x - 252-x + 420-x 


240 


2 1 

• x + 45 *x + 10) J 



4. Der Algorithmus von Casteljeau 

Mit diesem Algorithmus begannen die Arbeiten von de Casteljeau im Jahre 1959. Wir betrachten eine 
einfache Konstruktion zur Erzeugung einer Parabel. Die Verallgemeinerung dieser Konstruktion fuhrt 
zu den Bezierkurven und deren numerischer Berechnung. 

Es seien b 0 , b\ und b 2 drei beliebige Punkte des Anschauungsraums, die nicht auf einer gemeinsamen 
Geraden liegen. 

This algorithm was the start of de Casteljeau 's work in 1959. We look at an easy construc- 
tion of a parabola and its generalization will lead us to the construction of Bezier curves and 
their calculation. Let b 0 , bi and b 2 three points which are not lying on one common line. We 
set up the first two equations and substitute for b 0 , ^ and b 1j2 in the third one in order to obtain 
a parabola as a Bezier curve. 

This construction is a repeated linear interpolation. 0 < t < 1 gives the parabola within the 
triangle formed by b 0 , bi and b 2 . 

b ()i (t) = (\-t)b (l +tb ] 

b l2 (t) = ( 1 - t)b x +tb 2 teR 

^ 0 , 1,2 (0 — (1 — 0 \l 2 

b 0X1 (t) = (l-t) 2 b 0 + 2t(\-t)b x +t 2 b 2 , t^R. 

Wir setzen die beiden ersten Gleichungen in die dritte ein und erhalten mit dem quadratischen Term 
eine Parabel. 

Diese Konstruktion ist eine wiederholte lineare Interpolation. Wenn der Parameter t zwischen 0 und 1 
lauft, dann liegt die zugehorige Parabel innerhalb des von den drei Punkten b 0 , b\ und b 2 gebildeten 
Dreiecks. Dabei gilt boxiit = 0) = und 60 , 1 , 2 ^ = 1) = b 2 . 

Allgemein seien mit bo, b u b n insgesamt n+ 1 Punkte der Ebene gegeben. Wir bilden damit den 
Ausdruck 

K...i (0 = 0-0 K..J-1 (0 + 1 \ » 1 Zi^n. 

Dann ist 6o,i,...,«(0 der Punkt der Bezierkurve zum Parameterwert t g [0,1]. Das Polygon durch die 
Bezierpunkte bo, b\, b n i st das Bezierpolygon. Die Punkte der Bezierkurve zum Parameterwert t 
werden mit dem de Casteljeau Schema berechnet. 
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1 - 1 \ 


b b , b t , 

n n—l,n n-2,n-i,n 


-> *0.1 n=Pn( t \ t G [0,1] 


Bei der numerischen Berechnung wird nicht die vollstandige Dreiecksmatrix benotigt, denn die erste 
Spalte kann bei jedem Schritt geeignet uberschrieben werden. Es muss daher keine Tabelle abgespei- 
chert werden. 

(Zur Illustration kann der Algorithmus mit Schiilern handisch fur vier Punkte Z? 0 = (1,1), b\ = (1.5,0), 
b 2 = (3,1) und b 4 = (4,2.5) und t = 0,4 durchgefuhrt werden. Das lieBe sich natiirlich auch fiir den 
allgemeinen Parameter t durchfiihren. Josef) 

Given are n + 1 points b 0 , ..., b n . They form the Bezier polygon. The Casteljeau table 

returns the points of the Bezier curve for parameter t s [1,0]. This algorithm is more efficient. 
It is based upon a repeated linear interpolation. Neither equidistant nodes nor Bezier poly- 
nomials are needed. 

(As an illustration I'd recommend to perform this algorithm with pupils manually for eg four 
points: b 0 = (1,1), bi = (1 .5,0), b 2 = (3,1) and b 4 = (4,2.5) and t = 0.4. One could - should? - 
do this also producing the iterated curves. Josef) 
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Wir definieren die Funktion CB (p , t ) - sie liefert den Wert - und RB (p , t , a ) - sie liefert die a- te 
Spalte der Tabelle fur die Punktmatrix p. 

We define function CB(p,t) - returning the value - and RB(p,t,a) - returning the a th column 
of the table for given point matrix p. 

CB(p, t) := ITERATE (VECTOR ((1 - t)-x + t-x , j, DIM(x) -1), x, p, DIM(p) - 1) 

#23: j j + 1 

RB(p , t, a) := ITERATE (VECTOR ((1 - t)-x + t-x , j, DIM(x) -1), x, p, a) 

#24: j j + 1 


#26 : pi 


1 1 
1.5 0 

4 1 

#27: [RB(pl, 0.6, 0), RB(pl, 0.6, 1), RB(pl, 0.6, 2)] 

13 2 " 

10 5 IT 5S 13 

3 ' LL 25 ' 25 JJ 


#2S: 



" 1 

1 " 



3 

0 



2 



_ 4 

1 _ 



#29: [RB(pl, t, 0), RB(pl, t, 1), RB(pl, t, 2)] 


#30: 



" 1 1 " 


t 





+ 1 1 - t 



3 


2 

rr 2 2 ii 


0 

, 


, U_2-t + t + 1, 2-t - 2-t + lJJ 


2 


5 ■ t + 3 





t 



_ 4 1 _ 


2 



#31: CB(pl) 


r 2 2 1 

|_2 ■ t + t + 1, 2-t - 2 ■ t + lj 



#32: p2 := 


1 1 

1.5 0 

3 1 

4 2.5 

#33: [RB(p2 , 0.4, 0), RB(p2, 0.4, 1), RB(p2, 0.4, 2), RB(p2, 0.4, 3)] 

1 1 

r 1.2 0.6 " 

1.56 0.52 


#34: 


1.5 0 

3 1 

4 2.5 


" 1.2 

0.6 " 

2.1 

0.4 

_ 3.4 

1.6 _ 


2.62 0.SS 


, [[1.984, 0.664]] 


3 2 3 2 

3 ■ t - 6 ■ t - 3 ■ t - 2 3 ■ t - 12- t + 6 ■ t - 2 


#35: CB(p2 , t) = 


2 


2 
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#36: [RB(p2 , t, 0), RB(p2, t, 1), RB(p2, t, 2), RB(p2, t, 3)] 


#37: 



1 

1 



1.5 

0 



3 

1 

p 


4 

2.5 



0.5-t +1 1 - t 

1. 5 • (t +1) t 

t + 3 1.5-t + 1 


t + t + 1 


2-t - 2-t + 1 


- 0 . 5 • (t - 6-t - 3) 0 . 5 • t • (t + 4) 


IT 3 2 3 2 11 

, LL - 0. 5-(3-t - 6-t - 3-t - 2), - 0. 5-(3-t - 12-t + 6-t - 2)JJ 



Vergleiche mit den oben gefundenen handischen Ergebnissen - Compare with the results above 
derived manually 


5. Gradanhebung - Raising the degree 


Wir fiigen einen Punkt zum Bezierpolynom hinzu, ohne die Bezierkurve zu andem, wobei wir aber 
den Grad um 1 erhohen. Das ergibt neue Punkte die der Gleichung (1) geniigen. 


( 1 ) 




n+l 


2>, . a -ty- , t i ='Zb i l 

i = 0 V 1 J i=0 


n + l 


v * y 


(i — t) n+xi t‘ . 


Wir multiplizieren die linke Seite mit t + (1 - t) und erhalten durch Koeffizientenvergleich die Glei- 
chungen (2) und (3). 

^ + 1^ 


( 2 ) 

(3) 






b : + 


f n > 

yi ~ 1 j 


lr ! oder- or 


b}=- 


n + 1 


-b.+ 


1 — 


n + 1 


b, 


2 — 1 , 2 = 0 ,..., 72 + 1 


We add one point to the Bezier polygon without changing the Bezier curve. Doing this we 
raise the degree by 1. We find new points bi 1 which meet equation (1). Multiplying its left 
side by t + (1 - t) and comparing the coefficients leads to (2) and (3). The new points are 
emerging from the old ones by linear interpolation at the values i/(n+1). Raising the degree 
is important for numerical calculations because many programs require all Bezier polygons 
having the same degree (for exchangig data purposes). 
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Dann veranschaulichen wir an einem Beispiel die Polygone, die bei der Gradanhebung entstehen. 
Gegeben sei die Punktmatrix p3. Als neue Funktion wird MP (p) defmiert, die die Eckpunkte des 
neuen Polygons liefert. Dann zeichnen wir die zugehorige Parabel bzw. den Parabelteil und die Poly- 
gone, die bei der Gradanhebung entstehen. Diesen Vorgang konnen wir beliebig oft wiederholen 

MP(p) returns the vertices of the new polygon. Then we plot the parabola and the new poly- 
gon. This procedure can be repeated arbitrarily often. 




r r 

j 


3 "] 

■v 

3S: Ml 

B(p) := VE( 

ITOR IF j = 0 , 0, 


-p +i 

-p , o, i 

OIM(p) 



^ V 

DIM(p) 

1 i . 

o 

:m(p) J : + i 

j 


2 1 
0 4 

#40: GB(p3 , t) = [4-t-(l - t) , + 1)] 

" 0 0 

4 2 

3 3 

#41: (p4 := MB(p3)) = p4 := 

4 

2 

3 

_ 0 4 

#42: GB(p4 , t) = [4-t-(l - t) , 2-t-(t + 1)] 




Wir fahren fort bis zum Polygon berechnen wieder GB (pi 0 ) und stellen fest, dass sich immer 
wieder die gleiche Parabel als Bezierkurve ergibt. Alles gemeinsam kann in unterschiedlichen Fraben 
gezeichnet werden. 

We proceed until polygon plO, calculate again GB(plO) and observe that the result is the 
same parabola as Bezier curve. Finally we plot all together using different colours. 
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Alle Polygone, die jeweils durch Gradanhebung auseinander hervorgehen gemeinsam mit der Bezier- 
kurve - hier eine Parabel. 


6. Ausblicke - Prospects 

Es ist einfach, die Bezierkurven auf beliebige Intervalle [a,b] zu verallgemeinern. Ableitungen und 
Integrale von Bezierkurven konnen noch untersucht werden. Werden mehr als 10 Punkte zur Kon- 
struktion einer Bezierkurve verwendet, dann werden die Berechnungen zu aufwandig und es werden 
Splinekurven in Bezier-Form beniitzt. Dies und weitere Untersuchungen sind in dem Buch von Farin 
teilweise durchgefiihrt. 

It is easy to generalize Bezier curves for arbitrary intervals [a,b]. Derivatives and integrals 
can be investigated. If one uses more than 10 points for designing a curve then the calcula- 
tions will become too bulky and one has to use splines in Bezier form instead. These and 
other investigations can be found in Farin 's book. 

7. Literaturverzeichnis - References 

Farin G., Kurven und Fldchen im Computer Aided Geometric Design. Eine praktische Einfiihrung , 
Vieweg 1994 

Kose K., Schroder R., Wieiczek K., Numerik sehen und verstehen. Ein kombiniertes Lehr- und Ar- 
beitsbuch , Vieweg 1992 

More books (also treating splines): 

Futher W., Ohlsmann M., Mathematische Grundlagen der Computergraphik , Vieweg 1989 
Muller B. (editor), Techniken der Computergrafik , Franzis-Verlag 1994 
Schonefeld S., Numerical Analysis via DERIVE , Math Ware 1994 (only Splines) 
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Weber M., 3-D-Grafik Theorie und Praxis , IWT 1985 
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http : //www. geometry. at/ahsktn/bezier/praesentation_copyright .pdf 
http : // www . dma . uf g . ac . at/ app/link/ Grundlagen : 2D-Graf ik/ module/ 1307 9 ?step=3 
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In times of DERIVE 6.10 we can additionally introduce slider bars to visualize how the form 
of the Bezier polygon influences the form of the Bezier curve. You can compare with working 
with a graphics program. The origin is a fixed point of the polygon, the other three points can 
be varied by using the slider bars. Josef 


Working with slider bars for the coordinates: 

#50: [xl :=. yl :=. *2 :=. y2 :=, *3 :=, y3 :=] 

f 0 0 1 


#51: pts 


xl yl 
x2 y2 


L y3 J 


#52: GB(pts) 



See also DNLs # 53 and #66. 
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Anmerkungen zu Thomas Weths Kurven aus DNL#14 
Comments on Thomas Weth's curves in DNL#14 

Peter Baum, Kassel, Germany 


Peter Baum sent a contribution dealing with Thomas Weths'article in the DERIVE 
Newsletter #14. It is always a pleasure for the editor so see that the DNL is a 
source for the readers' own investigations. Peter sent the MTH-file with German 
comments and many TIF-files of the plots. I try to translate his comments connect- 
ing the DERIVE commands for our many non German speaking friends. Josef. 

The mapping described by Thomas Weth in DNL# 14 (inversion on a circle followed by an reflection 
wrt the origin), can easily be represented by means of polar coordinates. Even the mapping of circles 
is very clear. I restrict myself to the description, the reader might try to derive the formulae. I will 
start with the mapping of straight lines. 



Auxiliary point Q is obtained by inversion of P wrt to the circle with centre in the origin. 

The final mapping of P is F as result of a reflection with centre M. 

Let (s, a) the pre-image (= german: Urbild) of a point and (r, a) its image. The circle of inversion 
with radius R = a has its centre in the origin. Then we obtain the equation of the mapping described 
above 

2 

2-a 

#1: r = - s 

s 

If the picture of a point lies in the origin (pole) then r = 0, which leads to s = a^l. All these points 
form a circle concentric to the inversion circle with radius s = a-J 2. I will call this circle pole circle. 

#2: r = a 

#3: r = a-/2 

#2 and #3 give these two circles (in polar coordinates, of course). 
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1. The pre-image curve is a straight line parallel to y = 0 with distance t. 

Its equation in polar coordinates is given by #4 and simplifying #1 returns the equation of the 
image curve: 

t 

#4: s := 

COSO) 

2 

2 ■ a -COSO) t 

#5: r = 

t COSO) 

Thanks the slider bars we can now easily animate this mapping. 



The picture has a vertex if the line is a tangent of the pole circle. Then t = a yfl. By substituting 
for t = a\fl in #5 we obtain the equation of a kissoid. 


# 6 : 

#7: 

#S: 

#9: 


r = 72- a -COSO) - 

COSO) 

r = - J2 ■ a - SINO) ■TANO) 


r = -/2-a-COSO) 


^2- a 

COSO) 


Equation #7 is the well known form of a kissoid and by simplifying it - forcing the cosine - we 
can observe the identity of both forms. 
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2. Now the pre-image curve is a circle with its centre in M(f,0) and radius b. 

The equation of the circle in polar coordinates is given by expression #6: 

#1: LOADED: \DfD\DNL\DNL95\MTH19\D0T_CIRC.MTH) 

2 

2 ■ a 

#2: r = - s 

s 

#3: r = a 

#4: KSTRX(0, 0, a.^2, 30, 0.6) 

#5: r = a -42 

2 2 2 
#6: s := t.COS(ffl) + VO - t .SIN(a) ) 

#3 and #5 are inversion circle and pole circle. #4 is a function which plots the circle as 
a dashed line. This function is part of a little toolbox which will be described at the end 
of this contribution. Josef 
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Simplifying expression #2 using the appropriate Mode Settings will lead us to the general 
equation of the picture curve: 

#7: Trigonometry := Auto 




Substracting the right sides of #15 and #18 results in 0 verifying the identity. 

#15 (and #18 of course) is only valid for t^b- which is a circle passing the origin - The shape 
of the image curve depends on the position of the pre-image circle. It has a vertex if it is osculat- 
ing the pole circle. Then t = a42+b or t = aj2-b. In the first case we will find ’’tears” or 

’’drops” left of x = 0, in the other case ’’tears” right of x = 0 ( 2 b < aj 2 ) or cardioids (2b > as! 2) 
will appear. See the following plots. The radius of the inversion circle a = 2, the pole circle is 
the dashed line - the plot of #4 from above, the pre-image circle is in red and its image in blue. 
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For ploting the dashed circle one has to switch - for only this single plot to rectangular 
coordinates and set the parameter interval to [0,1]. All plots are done applying the 
slider bar feature of DERIVE 6. 



This is the plot with the sliders in any position. 

The plot performed with the DOS-versions didn't show any gaps. (I tried this with versions 
2. xx and 3.xx.) For circles which don't contain the origin within their circumference the recent 
Derive version shows gaps in the region of the osculating points between the circle and its 
tangents from the origin. The gaps are appearing in the mapping, too. 

(We will compare the representations given by other CASs at the end of this contribution.) 



t = 272 + 1.5; b = 1.5 



t = 2^2-l-,b = l 
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Later we will see that for 2b = aj2 the "tear" will degernate to a kissoid. For increasing b the 
kissoid will change to a cardiod. If b keeps increasing then the cardioid will shrink and for 
b = a \/2 degenerate to the origin, b > a V 2 will result in a cardiod right of x = 0. 



t = 2^2-2;b = 2 



t = 2^2 -2.5; b = 2.5 



i = 2V2-3.1;Z> = 3.1 


If the pre-image circles intersects the pole circle then the image curve has a double point in the 
origin. It has the shape of a leminscate or of a snail of Pascal. Because symmetric lemniscates 
are intersecting the x-axis symmetrically wrt the origin we substitute a = 0 and then a = n in 
expression #15 above in order to obtain #18 and #19: 
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#18: r = 


#19 : r = 


2 2 2 2 2 2 

(2-a - b + t ) ■ | b | - t-(2-a + b - t ) 

2 2 

b - t 

2 2 2 2 2 2 

(2-a - b + t }-|b| + t-(2-a + b - t ) 

2 2 

b - t 


#20: SOLVE 


2 2 2 2 2 2 
(2-a - b + t )'|b| - t-(2-a + b - t ) 

2 2 
b - t 


2 2 2 2 

#21: t = - 7(2 -a + b ) v t = -a + b ) v t = 0 


2 2 2 2 2 2 
(2 -a -b + t ) ■ | b | + t* (2 ■ a +b -t) 

, t 

2 2 
b - t 


The two values for r should be equal because of the given symmetry. Then we solve the equa- 
tion for the possible distances of the centre of the pre-image circle. 

t = 0 gives a circle concetric to the inversion circle. 



a = 2;b = 3.30; t = 0 


Substituting t = ^2 a 2 +b 2 in #15 and simplifying gives the equation of a lemniscate (#23) and 

if the tangents in the double point are perpendicular (b = a-J 2 ) then it will become the Lemnis- 
cate of Bernoulli. If b decreases then the lemniscate will become smaller. 

2 2 2 2 
#23 : r = - 2-^C(2.a + b 3-C0S(o) - 2-a ) 

2 

#24: r = - 2.^/2. |a| .^(2-COSCffl) -1) 



a = 2; b = 5; t = V33 
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a = 2;b = 2V2; t = 4 
Lemniscate of Bernoulli 



a = 2; b = 1; t = 3 


A kind of snail of Pascal appears if the pre-image circle intersects the pole circle and contains 
the origin. The general form of this curve is r = a + b cos a. It is a pity that it is not possible to 
obtain this curve by finding a special form of #18. 



a - 2; b = 3; t = 2 


3. If the pre-image circle passes the origin then t = b and the image curve contains infinity. 

For t = b the equation of this special family of circles is given by #27: 

2 2 2 
#25: ss := t-COS(ffl) + 7(b - t .SIN(tx) ) 

#26: ss := |b-C0S(o)| + b-COS(o) 
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#27: ss := 2.b.C0S(a) 

2 

2 ■ a 

#2S: r = ss 

SS 


2 

a 

#29: r = 2-b.C0S(a0 

b-COSO) 

Simplifying the equation of the mapping gives #29 which is the image of these special circles. 
It is remarkable that this equation is identical with expression #5 from page 35 - the image 
curves of straight lines - if b • t = - a 2 . That means: to each line exists one circle passing the 
origin with the same image curve. Geometrically you can obtain the centre of the circle by 
finding the pole of the line wrt the inversion circle and then reflecting it wrt the origin. The kis- 
soid, the strophoid and the trisectrix can be produced as pictures of circles. 




It is fascinating to watch the various changes from one type of curve to another by appropriate varia- 
tion of the parameters. You can hardly find limits for playfulness with DERIVE. ( The deficit of the 
weak representation of the polar forms is more than substituted by the advantage of intro- 
ducing slider bars, Josef in 2009 ) 

A final remark for teachers: surely plotting equations is a comfortable application of the computer. 
But primarly we have to train the students' ability to imagine the graph of a given equation in the 
various systems of coordinates and to recognize its dependency on the parameters of the equation. For 
this purpose you have to practise sketching the curves by hand, too. After having acquired some ex- 
perience the use of the computer is a meaningful enrichment and completion. Peter Baum. 
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Preparing this contribution for the DNL I reproduced Peter Baum's plots and I had the idea to 
represent the pole circle as a dashed line (for 0 < <P< 1 and in rectangular coordinates) and 
the pre-image circle as a bold line. Here are the two DERIVE - tools with three examples fol- 
lowed by the DERIVE screen. Josef 


#1: KDOP(m, n, r, sc 1) 


m + (r - 0 , 02 ■ sc} ■ COS (2 ■ 7T ■ tp) n + (r - 0 , 02 ■ sc) ■ SIN(2 -jr-cp) 


_ m + (r + 0 , 02 ■ sc} ■ COS (2 ■ 7T ■ tp) n + (r + 0 , 02 ■ sc) ■ SIN(2 -jr-cp) _ 


f 


r 2-77 ^ 


f 2,77 Y 


#2: KSTRXOn, n, r, t, p) := VECTOR 

m + r-COS 

(i + p-tp) 

, n + r ■ SIN 

(i + p-tp) 

■H 

1 

+-■ 

□ 



c t J 


L t J. 



#3: 

#4: 

#5: 

# 6 : 

#7: 


KD0P(2, 1, 1.5) 

KSTRX(-1, -1, 1.8, 40, 0.8) 
KSTRX(2.5, -1,5, 1.2, 20, 0.5) 
KSTRX(2 . 5 , -1.5, 1.23, 20, 0.5) 



Let's look how other Computer Algebra Systems of 2009 are performing the plots of the 
curves in polar representation: 

The Voyage 200 behaves pretty the same like DERIVE, we see the gaps. 
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TI-NspireCAS has no problems - like the V200 and the TI-92 to do the calculations but the 
plots are a bit worse because the endpoints of the gaps are connected - which has a 
mathematical explication but does not look fine. I introduced slider bars which are not avail- 
able on the TI-92 and Voyage 200: 



I chose the minimized form of the slider bars. 



That it can be done better shows the plot produced with WIRIS (including slider bars, too): 



s=tcos(a) + v b 2 -t 2 sin (a) 2 ; 

-2 t- cos (a) V-t 2 sin(a) 2 + b 2 + (-t 2 cos (a) 2 +t 2 sin (a) 2 ) + (2a 2 -b 2 ) 

r ™ * 7 

V _ t 2 sin(a) 2 + b 2 +t cos(a) 

a : =slider(O..5..0.2) ;b : =slider(0. .5. .0.1 ) ;t : =slider(— 10. .10. .0.5) ; 
plot({a,b,t}) — ► plotterl 
plot(polar_curve(a,a,0..2fi)) plotterl 
plot(polar_curve(a- /2 ,a,0..2fi)) plotterl 
plot(polar_curve(s,a,0..2?i) J {color=red}) plotterl 
plot(polar_curve(r,a ,0. .2n) ,{color=blue,line_width=3}) plotterl 
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Raymond Wieloszynski, Six-Fours les Plages, France 

Monsieur, 

I desire to pose a question, please excuse for my English: 

How to develop by series entieres in zero the next equation differentielle: 

4x (1 - x) y" + 2(1 - 3x) y' - y = 0 

or this equation differentielle in serie entiere solution of this equation y(0) = 1 

x y" + (1 - 2x) y' - 2y = 0 


I thank you for your reply. 


DNL: I am not an expert in this field but I looked into the manual of DERIVE 3.xx and found a function in the 
utility file ODEAPPR.MTH which might help you (even with DERIVE 2.50). 

I add a listing for the generalized problem and copy the respective explanation of the function from the DERIVE 
manual (now from the DERIVE 6.10 online help!): 

TAYLOR_ODE2 (r, x, y, v, xO, yO, vO, n) simplifies to an nth order truncated Taylor series solution of 
the equation y" = r(x, y, v). Here v represents y', and vO is the value of v at x = xO. For example, 
TAYLOR_ODE2 (2-x-v + x 2 y + 3 x, x, y, v, 0,0,1, 5) simplifies to 

5 3 

3 • x 5 • x 

+ + x 

10 6 

To verify that this is an approximate series solution to the differential equation, substitute this series 
solution for y and substitute the derivative of this expression for v into the unsimplified difference be- 
tween the two sides of the equation. Then use the Calculus Taylor Series command to verify that the 
third order Taylor series approximation of that difference is 0. The reason for using a third rather than 
fifth order Taylor series approximation is that the second derivative of the truncated series reduces its 
order by 2. 

TAYLOR_ODE2 simplifies to ? if r is not sufficiently differentiable for an nth order truncated series. 
Try a smaller value for n if this happens. 


Now let's have a try: 

f y - 2-(l - 3-x)-v 'l 

#1: TAYL0R_0DE2 , x, y, v, 0, yO, vO, 3 

k 4 • x • (1 - x) , 

#2: ? 

f y - 2- Cl - 3-x)-v 3 

#3: TAYL0R_0DE2 , x, y, v, xO, yO, vO, 3 

l 4-x-Cl - x) J 

3 2 2 2 

x -C2-v0-C14-x0 - 9-x0 + 3) + yO-C7-xO - 3)) - 3-x0-x • C2-v0- C20-x0 - 17-xO + 5) + yO-C9-xO - 5)) +~ 


2 2 2 2 
3-x0 •x-C2-vO-(34-xO - 41-xO + 15) + y0-(ll-x0 - 7)) - xO • (2-vO-xO- (56-xO - Sl-x0 + 33) - y0-(35~ 


2 2 
48-xO -CxO - 1) 

2 

•xO - 87-xO + 48)) 


# 4 : 
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The first question mark is the answer to use x = 0 but for x = 0 function r is not defined. So I did once more and 
generalized: y(xO) = yO and/(xO) = vO. 

Just to demonstrate how the function works I try with j/(x = 2) = 0, y'(x = 2) = 1. 1 ask for the series of order 5: 


y - 2- ci - 3 -x) -v 

#5: TAYL0R_0DE2| , x, y, v, 2, 0, 1, 5 

4-x-Cl - x) 


#6: 


5 4 3 2 

1535 -x - 17298-x + 79608 -x - 189136- x + 238128-x - 128928 


6144 


This was all what I could do in 1995. In the meanwhile Terence Etchells provided a function TaylorODE which 
can be found in the UsersYEquationSolving subdirectory of Derive: 

TaylorODE. dfw Defines the function TaylorODE(f, o, t, h, n, init vals, header?, xy only?) that solves initial 
value Ordinary Differential Equations of any order using Taylor’s method. 


#15: 


#16: 


yO - 2-0 - 3 -x) -yl 

Tay1orODE| , x, y, 2, 5, 0.1, 10, [2, 0, 1] 

4-x-Cl - x) 


jiitl. 


2 ] 


yO - 2- Cl - 3 -x) -yl 

Taylor0DE| , x, y, 2, 5, -0.1, 10, [2, 0, 1] | |u[l, 2] 

4-x-Cl - x) 


D, [2, 0, l]Jji 


I plot function #6 and the points created by #15 and #16 on the same axes: 



(The respective file is among the DNL#19-files.) 

Maybe that any of our members will find a more suitable answer? 
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ISAAC 1997 

From: rcorless@rawhide.cecm.sfu.ca (Rob Corless) 

Subject: For distribution to the DERIVE User Group 

We (Bruce Char and Rob Corless) are considering putting in a bid in to hold ISSAC 97 in Hawaii. There is 
little doubt in our minds that most potential attendees would find Hawaii a pleasant venue in which to hold a 
meeting. Our preliminary investigations into local arrangements and travel makes us believe it will be a satis- 
factory venue from those points of view as well. 

Our current thinking has centered on Kihei Maui, on the shore of southeastern Maui, a few miles from the 
Maui Supercomputing Center. We have a quote from the Maui Intercontinental Hotel of approximately $100 
US per night per room, single or double. There are accommodation alternatives (e.g. rental of condominiums 
that house four or six) in the area of Kihei as well. Living expenses are otherwise comparable to the continental 
USA. There are regularly scheduled commercial flights direct from the US mainland to Maui; ground transfer 
from the Maui airport to Kihei would be approximately 30 minutes via van. 

Airfare information at this early date is sketchy and vague; the current "21 day advance reservation" price is 
approximately $1100 US return from Philadelphia to Maui; from Vancouver, with a block booking of 100 
seats, the price is $499 (Canadian) return to Honolulu. Return airfare between Honolulu and Maui is approxi- 
mately $150 US. Qantas flies from Sydney to Honolulu for $1295 (Australian). The fact that Canadian prices 
appear to be the cheapest above may reflect the fact that a Canadian travel agent did most of the price search; 
cheaper flights may be available through your local agent. (For example last year during a sale, the return air- 
fare between Philadelphia and Honolulu was approximately $700). The agent was unable to give quotes for 
flights from Tokyo or from Europe, which would clearly be useful to know. 

Expressions of support have been received for our bid both by the Maui Supercomputing Center and the Soft 
Warehouse (creators of Derive). They would assist us in making local arrangements. The Supercomputing 
Center has extensive prior experience in setting up computing meetings at the Intercontinental Hotel. 

There is also the possibility of running another meeting (Stanly Steinberg's Applied Computer Algebra meet- 
ing) 'back to back' with ISSAC, to optimize the use of air transportation and clout with hotel bookings. We 
would appreciate suggestions for other meetings that we try to could co-schedule with ISSAC '97. 

We are also strongly interested in encouraging people who think of themselves as numerical analysts to 
come to this meeting. Some ISSAC meetings in the past have had a numerical flavour; with the Maui High 
Performance Computer Center supporting the meeting, it seems quite reasonable to expect that this meeting will 
have a similar sub-theme of numerical/symbolic interaction. 

Please take the time to respond to the following questionnaire. 

1. If ISSAC '97 was held on one of the Hawaiian islands, would you be interested in attending? 

2. Would you be able to come? 

3. Block-bookings or charter flights may significantly lower the cost of airfare. This requires planning ahead 
Would you be able to do so? 

4. The usual date for ISSAC is in the (northern hemisphere) summer. The climate of Hawaii is uniformly pleas- 

ant. and a summer conference there would provide a welcome change for those living in hot, muggy areas in 
(for example) the U.S. and Canada. However, several people have indicated informally that they would like 
to have an ISSAC meeting in Hawaii during the (northern hemisphere) winter, say January or February. How 
do you feel about that? 

5. Would you classify yourself as a numerical analyst or a computer algebraist (for want of a better name for 

people interested in symbolic computation) or both or neither? 

6. Please give us any other comments you have. 

Send all responses to bchar@mcs.drexel.edu or rcorless@uwo.ca 
Thank you for your attention 


Bruce Car and Rob Corless 


p46 


DERIVE POETRY 


D-N-L#19 


Dear Josef 

A student recently submitted an essay on the use of DERIVE in mathematics education. I have en- 
closed an amusing quote from his essay, that may appeal to readers of the DERIVE Newsletter. 

"HTtoteblebutn anil BEtoerbUbf e agreeb to (jabr a battle, 

Jfor CtoecbUtiee daib Cfaceblcbum tab spoileb tjid nice tieto rattle, 

Judt ttjen tletu bp a monstrous croto 

Sis ( big ad a tar barrel 

<C03t jut) trigbteneb bottj our tcroetf do 

lEtitp quite forgot i&eir quarrel. 4 * 

c^oii there In nutshell we have its acolytes 

depicted us a monstrous crow while the arguments for and 
against such primitive togs ns calculators- are cut off by the 
arrival of this terrifying creature. 

Best wishes 
Ted Graham 










